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Abstract 

We  have  prepared  nanoscopic  gold  cylinders  of  controlled  radius 
and  aspect  ratio  via  electrodeposition  of  the  metal  within  the  pores  of 
anodically-grown  porous  aluminum  oxide  membranes.  The  nanocylinder 
radii  are  determined  by  the  pore  dimensions  of  the  host  alumina  which, 
in  turn,  depend  on  anodization  conditions.  The  particle  aspect  ratios 
were  controlled  by  varying  the  amount  of  An  deposited  within  the  pores. 

The  optical  spectra  of  the  gold  nanocylinder/alumina  composites  exhibit 
strong  absorption  bands  in  the  visible  spectrum.  The  extinction  maxima 
(Xmax)  values  for  gold  particles  approaching  sphere-like  geometry  agree 
well  with  Mie  Theory  calculations.  The  blue  shift  of  Xmax  as  the  particle 
aspect  ratio  is  increased  is  in  qualitative  agreement  with  Maxwell-Gamett 
Theory.  We  propose  a  simple  modification  of  Maxwell-Gamett  Theory 
that  addresses  both  size  and  shape  effects. 


I.  Introduction 

The  optical  properties  of  nanoscopic  metal  particles  have  attracted 
a  great  deal  of  interest  for  many  decades.  Indeed,  it  has  been  almost  a 
century  since  Gustav  Mie  demonstrated  how  particle  size  gives  rise  to  the 
various  colors  of  colloidal  gold  solutions  [1.2].  In  more  recent  times 
attention  has  focused  on  the  application  of  such  small  metal  particles  in 
surface  enhanced  spectroscopy  [3],  photocatalysis  [4.5]  and  selective 
solar  absorbers  [6]. 

Most  work  to  date  on  the  optical  properties  of  nanoscopic  metal 
particles  has  involved  methods  that  preclude  control  of  particles  shape 


2 


and  orientation  with  respect  to  the  electric  field  of  the  incident  light. 
Exceptions  to  this  include  the  work  of  Liao  and  coworkers  who  prepared 
ordered  arrays  of  ellipsoidal  silver  particles  with  aspect  ratios  as  high  as 
3:1  using  photolithography  [7].  In  our  laboratories  we  have  developed 
methods  based  on  the  template  synthesis  strategy  in  which  various 
materials  are  deposited  within  the  pores  of  host  membranes  such  as 
track-etched  polycarbonate  [8]  and  anodic  aluminum  oxide  [9].  The 
microstructure  of  the  host  material  and  other  conditions  of  deposition 
thus  define  particle  geometry  and  orientation. 

In  our  efforts  to  understand  the  optical  properties  of  submicron 
metallic  particles,  we  have  found  that  anodic  aluminum  oxides  are  ideal 
host  templates:  they  are  optically  transparent  over  a  wide  range  (0.3  pm 
<  X  <  6  pm)  and  the  parallel  alignment  of  the  cylindrical  pores  means 
that  the  template-synthesized  metal  particles  are  also  aligned  in  one 
direction.  Hence,  the  particles  can  be  studied  without  removing  the  host 
template,  and  there  is  no  ambiguity  in  particle  orientation. 

In  our  early  work  we  relied  upon  commercially  available  anodic 
aluminas  [9].  which  contain  pores  ca.  0.25  pm  in  diameter.  In  these 
membranes  we  prepared  cylindrical  gold  particles  with  aspect  ratios 
(defined  as  the  particle  length.  A,  divided  by  its  diameter.  B)  as  high  as 
10:1.  The  resulting  Au/A1203  composites  showed  significant 

transparency  in  the  infrared  spectrum  but  not  at  shorter  wavelengths 
where  0.25  pm  diameter  particles  scattered  appreciably. 

In  a  recent  note,  we  presented  preliminary  results  of  our  optical 
characterization  of  gold  cylinders  template-synthesized  in  anodic 
alumina  membranes  prepared  in  our  laboratory  (10).  The  membranes 
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prepared  in  that  study  had  much  smaller  pores  (ca.  50  run  in  diameter), 
and  the  resulting  particles  showed  extinction  maxima  in  the  visible 
spectrum  reminiscent  of  aqueous  gold  colloids  of  similar  size.  In  that 
work,  we  demonstrated  how  the  template  synthesis  approach  allows  for 
simple  control  of  gold  nanocylinder  aspect  ratio  and  that  varying  the 
particle  aspect  ratio  gives  rise  to  systematic  shifts  in  the  extinction 
maximum  [10].  We  also  showed  that  the  observed  effects  of  particle 
shape  are  qualitatively  consistent  with  the  predictions  of  Maxwell- 
Gamett  Theory  (M-G),  which  is  a  simple  treatment  for  composite 
materials  that  addresses  particle  shape  and  orientation  but  not  size  [11]. 

In  this  paper  we  present  the  results  of  a  study  in  which  we  have 
systematically  varied  both  nanocylinder  aspect  ratio  and  diameter.  We 
show  how  the  template  synthesis  method  is  a  convenient  means  by 
which  one  can  control  these  parameters  and  hence,  the  composite's 
optical  properties.  Because  there  is  no  ambiguity  in  nanocylinder 
orientation  with  respect  to  the  incident  field,  the  analysis  of  size  and 
shape  effects  on  optical  spectra  should  be  greatly  simplified.  We 
compare  the  spectra  of  the  gold  nanocylinder/ alumina  composites  with 
the  predictions  of  M-G  Theory  and  Mie  Scattering  Theory  [2].  We  also 
propose  an  approximation  within  the  M-G  framework  that  addresses 
both  the  size  and  shape  of  the  nanoscopic  particles. 
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II.  Theory:  Mie  Scattering  and  lfaxwell-Garnett  Treatments  for 

Nanoscopic  Particles 

In  this  section  we  review  two  approaches  that  we  use  to  model  the 
optical  spectra  of  composites  containing  nanoscopic  particles.  Mie 
Scattering  Theory  is  a  rigorous  formulation  for  particles  of  any  size  but  is 
limited  to  spherical  geometry  and  interparticle  separation  distances 
which  are  much  larger  than  the  wavelength  of  the  incident  light  [1,2]. 
Maxwell-Gamett  (M-G)  Theory  is  strictly  valid  only  in  the  limit  where 
particle  dimensions  and  separation  distances  are  infinitely  small  relative 
to  the  wavelength  but  can  be  easily  generalized  to  particles  of  various 
shapes  [11].  In  our  experiments,  because  we  vary  both  size  and  shape  of 
the  gold  particles,  each  of  these  approaches  has  something  to  offer  in  the 
interpretation  of  the  optical  spectra. 

In  the  Mie  approach,  the  total  transmittance  through  a  film 
containing  spherical  metal  particles  is  given  by  [2a] 

Ttot =  exp  ( -  N  ?t  a  Qextd  )  (1) 

where  N  is  the  number  concentration  of  spheres  per  unit  volume,  a  is  the 
sphere  radius,  and  d  is  the  film  thickness.  The  extinction  coefficient  Qext 
is  related  to  the  Mie  scattering  coefficients  an  and  bn  through  [2a] 

oo 

Qext  =  fell  ( 2n+l )  Re  ( a„+bn )  (2) 

n=l 
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where  x  =  2nano  A  (no  is  the  refractive  index  of  the  host  medium  and  X  is 
the  wavelength  of  the  incident  light  in  vacuo).  The  term  Re  in  Equation  2 
signifies  the  real  part  of  the  complex  function  in  parentheses.  The 
scattering  coefficients  an  and  bn  describe  the  amplitudes  of  the  scattered 
electric  and  magnetic  fields,  respectively  [10],  and  are  given  by 

=  V.(y)  ¥,00  -  my„(y)V,.«  (3a) 

¥Vy)C«00  -  |"W)  C',00 


^  =  myn(y)Vi,(x)  -  ¥n(y)¥,M  (3b) 

m  V'nCy)  Cato  -  ¥,&)  C',00 

where  m  is  related  to  the  complex  refractive  index  of  the  metal.  nm,  and 
to  the  host  medium  refractive  index,  no.  via 

m  =  Sm/iio  (4a) 

and 

•  Do  =  Dm  +  i  km  (4b) 

where  nm  is  the  real  refractive  index  of  the  metal  and  km  is  the 
absorption  coefficient.  The  quantity  y  in  Equations  3a  and  3b  is  equal  to 
The  primes  In  Equations  3a  and  3b  signify  first  derivatives 
with  respect  to  the  argument. 

The  functions  \|fn  and  £n  are  related  to  the  Bessel  (JJ  and  second 
order  Hankel  (1$)  functions  via  [2, 12]. 
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y„(z)  =  (z  %l2)xn  j-w(z) 


(5a) 


and 

Cn(z)  =  (zJC/2),/2H&,(z)  (5b) 


where  z  is  equal  to  x  or  y. 

For  metal  particles  less  than  col  100  run  in  diameter,  the  Mie 
scattering  treatment  for  visible  wavelengths  requires  only  the  first  two  or 
three  expansion  terms  in  Equation  2.  Nonetheless,  the  computation  is 
still  quite  tedious.  A  quite  different  approach  which  is  valid  when  the 
metal  particle  dimensions  are  much  smaller  than  the  wavelength  of  light 
was  developed  by  J.C.  Maxwell-Gamett  [9a].  M-G  Theory  assumes  that  it 
is  reasonable  to  define  an  effective  complex  dielectric  constant  £c  for  a 
composite  containing  metal  particles  (with  complex  dielectric  constant 
Em)  imbedded  in  a  host  matrix  (with  complex  dielectric  constant  Eo) 
through 


§£  ~  So  = 

\£c  +  2£o 

ten,  +  2eJ 

(6) 


where  fm  is  the  volume  fraction  of  metal  in  the  composite.  The  complex 
dielectric  constant  is  related  to  the  optical  refractive  index  n  and 
absorption  coefficient  k  via 
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e  =  n2  =  e'  +  i  e" 

(7a) 

E’  =  n2  -  k2 

(7b) 

E"  =  2nk 

(7c) 

where  e*  is  the  real  part  of  the  dielectric  constant  and  E"  is  the  imaginary 
component  (the  refractive  index  n  is  not  to  be  confused  with  the 
summation  index  in  Equations  2  and  3). 

Maxwell- Garnett's  original  work  addressed  spherical  inclusions, 
and  it  is  for  this  specific  geometry  that  the  factor  2  appears  in  the 
denominators  of  Equation  6.  This  factor  is  also  called  the  screening 
parameter  [1  lb]  and  is  in  general  denoted  K.  Hence,  while  K  =  2  for 
spheres,  long  needle-like  particles  oriented  with  their  axes  of  revolution 
parallel  to  the  direction  of  light  incidence  have  a  k  value  that  approaches 
unity.  On  the  other  hand,  K  approaches  infinity  for  thin  flat  discs 
oriented  with  their  axes  of  revolution  perpendicular  to  the  light  incidence. 

From  a  knowledge  of  Em,  Eo.  fm  and  K,  one  can  derive  the  effective 
dielectric  constant  £c  for  the  composite,  and  then,  via  Equations  7a-c,  the 
effective  composite  optical  constants  xic  and  kc  can  be  calculated.  The 
total  transmittance  is  then  given  by  113] 


Tiot  — 


1  - 


(nc- 1  )2  +  k? 


(nc+1)2  +  k?J 


exp  (-4rckcd/X)  (8) 


where  d  is  the  film  thickness  and  X  is  the  incident  wavelength.  The  term 
preceding  the  exponential  accounts  for  reflective  losses  at  the  film/air 
interface  [14]. 
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m.  Experimental  Section 

A.  Porous  Alumina  Membrane  Preparation 

Nanoporous  alumina  template  membranes  were  prepared  by 
anodization  of  aluminum  (Aldrich  2mm  plate  99.999%  purity)  in  aqueous 
oxalic  acid  media.  Anodization  of  aluminum  to  yield  porous  oxides  is  a 
well  established  technique  and  has  been  discussed  extensively  in  the 
literature  [15].  For  details  on  how  the  porous  oxide  is  removed  from  the 
aluminum  substrate,  consult  Reference  16. 

Our  two  electrode  anodization  cell  consisted  of  an  insulated 
thermostat-controlled  plastic  bucket  (ca.  19  L),  a  6  cm  x  10  cm  stainless 
steel  cathode  and  a  10  cm  x  10  cm  A1  plate  anode.  Voltage  control  was 
provided  by  a  Sorensen  DCR-9B  power  supply. 

Pore  diameters  in  anodic  aluminas  are  roughly  proportional  to  the 
applied  voltage  in  potential-controlled  anodization  [161.  In  this  study,  we 
prepared  membrar  _a  pores  of  four  different  diameters.  Anodization 
at  40  V  and  50  V  in  4%  (w/v)  oxalic  acid  yielded  membranes  with  pore 
diameters  of  approximately  60  nm  and  75  nm  respectively.  Anodization 
at  70  V  and  90  V  in  2%  (w/v)  oxalic  acid  produced  pores  with  diameters 
of  approximately  86  nm  and  120  nm  respectively.  The  electrolyte 
temperatures  was  maintained  between  0°  C  and  7°  C.  Membranes  were 
grown  from  30  pm  to  40  pm  in  thickness  before  removal.  Figure  1  shows 
transmission  electron  micrographs  (TEMs)  of  facial  sections  of  alumina 

k 

composite  membranes  prepared  under  the  above  conditions.  The 
porosity  of  these  membranes  is  in  the  range  of  25%  to  30%. 
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B.  Gold  Nanoparticle  Preparation 

Gold  nanoparticles  were  prepared  by  electrodeposition  of  Au(I) 
(Technic  Inc.,  Orotemp  24  ™  gold  plating  solution)  within  the  pores  of  the 
alumina  membranes.  The  details  of  this  procedure  are  given  in  Reference 
9b.  Four  amounts  of  gold  (0.20  C,  0.40  C,  0.70  C  and  1.00  C),  in  terms 
of  coulombs,  were  deposited  in  each  of  the  four  different  pore  diameter 
membranes  over  a  3.3  -  3.5  cm2  area,  depending  on  the  measured 
porosity.  The  potentiostat  employed  was  an  EG&G  Princeton  Applied 
Research  Model  1 73  equipped  with  a  Model  1 79  digital  coulometer. 

C.  Anodic  Alumina  and  Gold  Nanqparticle/Alumina  Composite 

Characterization 

Alumina  pore  radii  and  gold  nanoparticle  dimensions  were 
determined  by  means  of  transmission  electron  microscopy  ( JEOL  2000 
TEM).  The  thin  sample  sections  required  for  these  analyses  were 
obtained  by  diamond-knife  microtome.  UV/Visible/Near-IR  spectra  were 
obtained  using  a  Hitachi  U-3501  spectrometer.  Figure  2  shows  a 
schematic  of  the  experimental  optical  configuration  used  to  examine  the 
gold  nanocylinder/alumina  composites. 

IV.  Results 

The  template  synthesis  method  allows  for  straightforward  control 
of  nanocylinder  size,  shape  and  orientation.  Figure  3  shows  TEM  images 
of  cross  sections  of  alumina  membranes  (prepared  at  40  V  in  4%  oxalic 
acid)  after  the  deposition  of  Au  into  the  pores.  At  low  deposition  levels 
like  0.020  C  (0.06  C  Au/cm2),  the  particles  are  short  and  sphere-like.  At 
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higher  deposition  levels  such  as  1.00  C  (0.30  C  Au/cm2),  the  particles 
become  longer  and  needle-like  within  the  cylindrical  radius  held  constant 
by  the  pore  walls.  Table  I  summarizes  the  results  of  measurements  of 
particle  dimensions  for  the  various  pore  radius/gold  deposition 
combinations  as  obtained  from  TEM  images. 

The  TEM  images  of  the  gold /alumina  composites  also  show  that 
the  gold  particles  do  not  lie  at  a  uniform  distance  from  the  membrane 
face  but  rather  through  a  layer  ca.  1  pm  in  thickness.  We  discussed  the 
reasons  for  this  in  a  previous  paper,  and  also  how  this  non-uniform 
particle  placement  affects  the  optical  path  length,  d.  and  metal  volume 
fraction,  fm  [9b].  We  will  discuss  this  point  further  in  the  next  section. 

Figure  4  shows  UV/Vis/NTR  extinction  spectra  [17]  for  three 
composites  of  different  pore  radii  in  which  the  gold  nanoparticles  are 
short  and  nearly  spherical  (Curves  1,2  and  3).  In  each  case  the 
composite  exhibits  a  significant  transparency  in  the  NIR  region  of  the 
spectrum  and  a  strong  extinction  band  in  the  visible.  The  wavelength  of 
maximum  extinction  increases  with  increasing  particle  radius  giving  rise 
to  color  variations  in  the  different  gold  particle  radius  composites.  The 
transmitted  light  for  the  30  nm  radius  particles  appears  red-purple  while 
60  nm  radius  particles  appear  blue-green.  This  trend  is  consistent  with 
the  colors  exhibited  by  colloidal  gold  particles  of  different  sizes  in 
aqueous  solutions^].  The  60  nm  particle  composite  also  shows  a  weak 
shoulder  (at  ca.  550  nm)  on  the  broad  band  centered  at  655  nm. 

The  dependence  of  the  optical  spectra  on  particle  shape  is  shown 
in  Figure  5.  As  the  aspect  ratio  increases,  the  extinction  intensity 
increases  and  Xmax  shifts  to  shorter  wavelengths.  This  trend  is  seen  for 
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all  composite  systems  considered  in  this  study  and  is  summarized  in 
Figure  6.  where  L'.»  is  plotted  against  the  natural  logarithm  of  the 
aspect  ratio,  ln(A/B).  The  two  smallest  pore  radius  systems  (30  nm  and 
38  nm)  show  values  that  converge  at  high  aspect  ratios.  The  larger 
pore  systems  exhibit  similar  trends  with  respect  to  shape,  but  the 
values  are  red-shifted  relative  to  the  small  particle  radius  composites. 

Thus  there  are  two  trends  immediately  apparent  for  gold 
nanocylinder  particles  aligned  with  their  principal  axes  parallel  to  the 
direction  of  light  incidence.  First  of  all.  red-shifts  with  increasing 
cylinder  radius  for  particles  of  all  aspect  ratios  (Figure  6).  Secondly,  for 
constant  cylinder  radius.  Xmax  blue-shifts  as  the  particle  aspect  ratio  is 
increased  (Figure  6). 

V.  Discussion 

A.  Comparison  of  Spectra  qf  Sphere-Like  Gold  Particles  with  Me 
Calculations 

As  it  was  mentioned  in  the  previous  section,  deposition  of  small 
amounts  of  gold  gives  rise  to  particles  whose  aspect  ratios  are  nearly 
unity.  It  is  of  interest  then  to  compare  the  experimental  spectra  for  these 
spherical  particle  composites  with  spectra  calculated  using  Mie  Theory. 
The  calculation  requires  a  knowledge  of  the  particle  radii  (which  we 
obtain  from  the  TEM  measurements),  the  values  of  the  optical  constants 
for  the  metal  particles  and  host  membrane  for  the  spectral  range  of 
interest,  and  finally  the  particle  density  N  and  optical  path  length  d.  The 
values  of  nm  and  km  for  gold  were  taken  from  Reference  18.  For  the  host 
medium,  we  assumed  no  =  1.6  [19],  While  the  anodic  aluminas  used  in 
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this  study  arc  slightly  absorbing  in  the  visible  spectrum  (ko  <  10'3),  we 

assumed  ko  =  0  to  avoid  complex  arguments  in  the  parameter  x  of 

«■ 

equations  3a  and  3b  and  thus  simplify  the  calculation. 

We  computed  Qat  using  Equations  2-4  with  n  =  1  through  3.  We 
have  found  that  the  third  term  in  Equation  2  is  negligible  for  all  particle 
sizes  considered  in  this  study.  The  total  transmittance  was  then 
calculated  from  Equation  1.  with  particle  density  N  and  composite  path 
length  d  chosen  so  as  to  yield  extinction  intensities  comparable  to 
experiment  ( vide  infra). 

The  Mie-calculated  spectra  shown  in  Figure  7  are  in  good 
qualitative  agreement  with  the  experimental  spectra.  As  the  sphere 
radius  increases,  the  extinction  maximum  shifts  to  longer  wavelengths 
and  decreases  in  intensity.  The  calculated  spectra  also  exhibit  increased 
band  broadening  with  particle  size,  though  not  to  the  extent  seen  in 
experiment  (see  Figure  4). 

As  shown  in  Figure  6.  the  lowest  aspect  ratio  particles  are  still 
somewhat  prolate  (ie.,  their  aspect  ratios  are  greater  than  unity).  Hence, 
in  order  to  compare  the  experimental  and  Mie  Theory  values,  we 
extrapolated  the  experimental  L«  values  shown  in  Figure  6  to  unit 
aspect  ratio  using  a  quadratic  fit  of  Xmas  versus  ln(A/B).  Table  n 
summarizes  these  Xmax  values  extrapolated  from  experiment  as  well  as 
those  calculated  from  Mie  Theory.  It  is  evident  that  the  quantitative 
agreement  is  quite  satisfactory. 

The  Mie-calculated  spectrum  for  the  60  nm  radius  particles 

¥ 

exhibits  a  shoulder  at  530  nm  which  compares  favorably  with  the  weak 
shoulder  seen  in  the  experimental  spectrum  at  ccl  550  nm.  In  the  Mie 
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treatment  this  band  is  due  to  the  gold  sphere  quadrupole  mode  and 
arises  from  the  n  =  2  index  term  in  Equations  2  and  3  [20].  This  mode  is 
negligible  in  all  but  the  60  nm  radius  calculation.  For  the  particle  radius 
range  considered  in  this  study,  the  leading  dipole  term  (n  =  1) 
predominates. 

In  principle,  it  should  be  possible  to  determine  the  particular 
concentration  N  using 

N  =  3  fm  /  4  7t  a3  (9) 

As  we  discussed  in  a  previous  paper,  because  of  the  geometry  of  the 
oxide  pores  and  the  manner  in  ■which  the  particles  are  created  in  the 
pores,  the  metal  volume  fraction  should  be  related  to  the  intrinsic 
porosity  P  of  the  host  membrane  and  the  optical  path  length  d  via  [9b] 

fm  =  P(2  a/d)  (10) 

In  fact,  when  we  attempted  spectral  calculations  using  N  based  on 
Equations  9  and  10,  with  estimates  of  P  and  d  from  TEM  images,  the 
extinction  intensities  were  much  higher  than  those  seen  experimentally 
(the  Xmax  results  are  not  affected  by  N  in  the  Mie  treatment).  Although 
the  reasons  for  this  apparent  discrepancy  are  not  explicitly  known,  the 
values  of  the  experimental  and  simulated  spectra  converge  if  lower 
porosities  are  input  into  the  Mie  calculations. 
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B.  Particle  Shape  effects  in  the  Maxwell-Gamett  Limit 

If  a  composite  contains  particles  whose  dimensions  are  much 
smaller  than  the  wavelength  of  the  incident  light,  then  the  Maxwell- 
Gamett  formulation  of  Equation  6  is  applicable.  Because  Equation  6 
does  not  address  particle  size,  it  cannot  be  used  to  model  size  effects  on 
extinction  maxima  as  in  the  case  of  Mie  Theory.  On  the  other  hand, 
particle  shape  determines  the  screening  parameter  K,  and  hence,  will  be 
manifest  in  the  composite  optical  properties  ric  and  k«. 

In  our  attempts  to  model  shape  effects  in  the  experimental  spectra, 
we  treat  the  gold  nanocylinders  as  ellipsoids  of  revolution.  The  long  axis, 
A,  of  the  particle  is  thus  identified  with  twice  the  semimajor  axis  a  of  an 
ellipsoid.  The  semiminor  axis  b  is  then  the  particle's  cylindrical  radius 
where  2b  =  B,  the  short  axis  or  diameter  of  the  particle. 

The  M-G  screening  parameter  K  is  related  to  the  Lorentz 
depolarization  factor  q  through  [1  lb] 

K  =  I  -  1  (11) 

q 


The  factor  q  depends  on  the  shape  of  the  particle.  For  an  electric  field 
incident  along  the  long  axis  A  of  a  prolate  ellipsoid,  the  depolarization 
factor  q*  is  given  by  (2] 


where 


Qa 


1 


(12a) 
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C-  1  - 


(12b) 


In  our  study,  the  gold  particles  are  aligned  with  their  principal  axes 
parallel  to  the  direction  of  light  incidence.  The  electric  field  is  thus 
incident  along  the  short  axes  of  the  particles.  For  the  special  case  of  an 
ellipsoid  of  revolution,  we  can  relate  the  depolarization  factor  along  the 
long  axis  (q*)  to  qj>,  the  depolarization  factor  along  the  short  axis  via  [21] 


qb  = 


i-q« 

2 


(13) 


In  our  experiments,  low  amounts  of  deposited  gold  leads  to 
particles  that  are  nearly  spherical.  In  the  spherical  limit  (as  A/B-*l),  the 
depolarization  factor  qj>  is  equal  to  1/3  and  K  is  equal  to  2.  However,  as 
we  increase  the  amount  of  gold  deposited,  the  aspect  ratio  increases,  and 
as  it  approaches  infinity  (A/B  -»  «,).  Equations  1 1-13  lead  to  qt,  *  1/2 
and  K=l.  It  is  thus  relevant  to  examine  the  M-G-calculated  spectra  for 
1  <  K  £  2. 

Figure  8  shows  absorption  spectra  calculated  using  M-G  Theory  for 
composites  containing  gold  ellipsoids  of  varying  aspect  ratios  such  that 
1  <  K  £  2.  As  in  the  case  of  the  Mie  calculations,  the  optical  properties  of 
gold  were  taken  from  Reference  18.  For  the  host  oxide,  we  used  complex 
Do  values  obtained  in  this  laboratory  [22].  We  modeled  the  experimental 
membranes  as  three-layer  systems  in  which  the  gold/alumina  composite 
layer  is  sandwiched  between  two  layers  of  alumina  which  contain  no  gold 
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(see  inset  of  Figure  8).  Thus  there  are  four  interfaces  where  reflective 
losses  occur  in  addition  to  the  absorption  losses  in  the  composite  and 
alumina  layers.  (9b]. 

The  M-G-calculated  spectra  shown  in  Figure  8  agree  qualitatively 
with  experimental  data  (see  Figures  5a  and  5b).  As  the  gold  particle 
aspect  ratio  increases,  the  absorption  intensity  increases  (due  to  the 
attendant  increase  in  path  length  and  metal  volume  fraction).  However, 
the  increase  in  aspect  ratio  also  engenders  a  decrease  in  K  and  Xmax 
shifts  to  shorter  wavelengths.  This  effect  is  seen  in  the  experimental 
spectra  for  all  pore  diameter  systems,  as  summarized  in  Figure  6.  The 
lowest  curve  in  Figure  6  describes  the  shape  effect  on  W  in  the  M-G 
limit.  All  experimental  Xmax  values  are  shifted  red  relative  to  the  M-G 
limit. 

In  the  M-G  limit  the  dependence  of  Xmas  on  particle  shape  can  be 
understood  in  terms  of  the  particle  dipole  resonance  condition  often 
discussed  in  the  surface  enhanced  Raman  spectroscopy  (SERS) 
literature  [3].  The  wavelength  of  highest  absorption  in  the  metal  particle 
composite  corresponds  to  the  condition  where  the  induced  field  in  the 
metal  is  maximal.  This  occurs  when  the  real  component  of  the  metal 
dielectric  constant  E'm  is  equal  to  -K£'0.  Hence,  the  right  hand  side  of 
Equation  6  becomes  very  large  (assuming  the  imaginary  component  of 
the  metal  and  host  dielectric  functions  are  very  small  and  the  composite 
dielectric  function  is  maximized).  The  real  dielectric  constant  £'m  for  gold 
(and  most  metals)  is  negative  in  sign  and  monotonically  decreases  at 
wavelengths  above  the  plasma  edge  [23].  Hence,  if  K  decreases,  the 
dipole  resonance  condition  occurs  at  a  shorter  wavelength. 
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C.  An  Approximate  Dynamical  McucweU-Qamett  Model  for 
Addressing  Size  and  Shape  Effects 

Mlc  Theory  successfully  describes  size  effects  but  is  limited  to 
spheres.  In  contrast  M-G  Theory  predicts  qualitative  trends  in  Xmi.*  with 
shape,  but  cannot  deal  with  size.  While  there  have  been  numerous 
scattering  treatments  for  ellipsoidal  [24]  and  cylindrical  particles  [25], 
they  are.  unfortunately,  even  more  tedious  than  the  Mie  treatment  for 
•spheres.  It  is  desirable  to  have  a  simple,  if  very  approximate,  model 
which  deals  with  size  and  shape  effects  simultaneously,  and  is  also 
appropriate  for  systems  in  which  interparticle  separations  cannot  be 
considered  infinite,  as  is  formally  required  in  single  particle  scattering 
treatments  such  as  Mie  Theory.  In  this  section  we  propose  such  a  model. 
The  basic  strategy  is  to  incorporate  a  size-  d  ependent  depolarization 
factor  into  the  M-G  treatment  and  then  calculate  an  effective  composite 
dielectric  function  from  which  we  can  simulate  spectra. 

In  1983,  Meier  and  Wokaun,  in  an  effort  to  predict  metal  particle 
size  effects  in  SERS,  derived  an  alternate  expression  for  the 
depolarization  factor.  Specific  for  spheres,  their  model  defines  an 
effective  depolarization  factor  q,fy.  which  is  related  to  sphere  radius  s  via 
[26] 


qeff  =  J  -  (  j)  k2  a2  -  (!j  i  k3  a3  (14) 

where  the  first  term  on  the  r.h.s.  of  Equation  14  is  the  Lorentz 
depolarizing  factor  (q  =  1/3)  for  spheres.  The  second  term  involves 
dynamic  depolarization  and  accounts  for  the  fact  that  a  particle  that  is 
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not  small  relative  to  the  Incident  wavelength  will  experience  different 
phases  of  the  Incident  field  [27].  The  term  k  is  equal  to  2ia\0/’k  (and 
should  not  be  confused  with  the  absorption  coefficient  of  Equations  4.7 
and  8).  The  third  term  accounts  for  damping  of  the  induced  dipole  due 
to  radiation  emission  [27]. 

The  expression  for  spheres  in  Equation  14  was  generalized  to 
ellipsoids  of  revolution  by  Zcman  and  Schatz  (27,28).  The  dimension  of 
relevance  in  the  size  dependent  terms  is  the  ellipsoid  axis  parallel  to  the 
electric  field  [26-28].  Therefore,  we  can  calculate  an  effective  screening 
factor  Keff  using 

Q«rr  =  qb  -  (j)  k2  b2  -  (|)  i  k3  b3  (15a) 

Kcff  =  —  -  1  (15b) 

qeff 

where  qj>  is  the  depolarization  factor  for  infinitely  small  particles,  and  b  is 
the  semiminor  axis  of  the  particles.  The  choice  of  b  follows  from  the 
orientation  of  our  particles  in  the  optical  experiments  (see  Figure  2). 

Using  the  same  approach  as  in  Section  VB  but  with  K^r  replacing  K 
in  Equation  6,  we  calculated  spectra  for  different  nanoparticle  radii  and 
aspect  ratios.  Figures  9A  and  9B  show  these  calculated  spectra  for  the 
30  run  and  60  nm  particle  radius  cases  respectively-  In  each  case,  as  the 
particle  aspect  ratio  increases,  shifts  to  shorter  wavelengths.  Also, 
as  seen  clearly  in  the  experimental  results,  the  extinction  bands  for  the 
60  nm  set  is  red-shifted  relative  to  the  30  nm  radius  case,  and  the  bands 
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are  significantly  broadened  (again,  in  accord  with  experiment).  Figure  10 
summarizes  the  calculated  values  as  a  function  of  aspect  ratio.  The 
general  similarity  to  Figure  6  is  quite  apparent  particularly  with  regard 
to  the  degree  to  which  the  different  pore  radius  systems  are  sensitive  to 
shape.  For  example,  the  experimental  Xmmx  values  for  the  smallest 
particle  systems  (30  nm  radius)  shift  about  20  nm  going  from  spherical 
(A/B  »  l)  to  needle-like  (A/B  -» «*)  geometry.  On  the  other  hand,  the 
largest  systems  (60  nm  radius)  blue-shift  ca.  70  nm  going  from  spherical 
to  needle-like  shapes.  The  analogous  shifts  based  on  the  dynamical  M-G 
calculation  are  48  nm  for  the  30  nm  radius  set  and  132  nm  for  the 
60  nm  systems. 

The  dynamical  M-G  approach  overestimates  L«  in  all  cases  and 
most  seriously  at  the  largest  particle  radius  (for  the  spherical  case,  see 
Table  II).  This  is  not  terribly  surprising,  since  we  know  from  the  Mie 
calculations  that  60  nm  radius  spheres  scatter  considerably.  Indeed, 
what  we  are  attempting  to  do  in  the  dynamical  M-G  model  is  to  fold 
scattering  losses  into  the  composite  absorption  coefficient,  kc.  For  three 
of  the  particle  size  regimes  examined  in  this  study,  this  approach  is 
certainly  not  quantitatively  correct.  However,  it  does  offer  simplicity  and 
a  reasonable  description  of  trends  in  band  shape  and  band 
intensity. 

VI.  Conclusions 

We  have  used  aluminum  anodization  and  the  template  synthesis 
method  to  prepare  nanoscopic  gold  particles  of  various  radii  and  aspect 
ratios.  Optical  spectra  for  the  gold/alumina  composites  containing 
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spherical  gold  particles  are  quantitatively  amenable  to  Mie  Theory.  The 
spectra  of  the  higher  aspect  ratio  particle  composites  show  both  a  size 
and  shape  dependence,  the  latter  following  qualitatively  the  predictions 
of  Maxwell-Gamett  Theory.  We  propose  a  simple  modification  of  the  M-G 
treatment  using  dynamical  models  for  the  depolarization  factor  q 
developed  for  SERS  theory  {26-28]  to  deal  with  size  and  shape  effects  in  a 
unified  calculation.  While  this  dynamical  M-G  approach  overestimates 
Xmax.  it  does  capture  qualitative  trends  in  Xmax.  as  well  as  band  shape 
and  intensity  as  the  size  and  shape  of  the  metal  particles  is  varied. 
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Table  I. 

Experimental  Results:  Gold  Nanoparticle  Dimensions 


Anodic 

Voltage 

(Volts) 

No.  C 
And)* 
(Coulombs) 

Particle 

Radius* 

(nm) 

Particle 

Length* 

(nm) 

Aspect 

Ratio6 

(A/B) 

K* 

40 

0.20 

30  ±3 

92  ±  12 

1.53 

1.65 

40 

0.40 

•i 

132  ±6 

2.20 

1.46 

40 

0.70 

H 

2681 16 

4.47 

1.22 

40 

1.00 

»f 

434  1  22 

7.23 

1.14 

50 

0.20 

38  ±2 

9818 

1.29 

1.78 

50 

0.40 

•« 

16016 

2.11 

1.48 

50 

0.70- 

if 

25617 

3.37 

1.30 

50 

1.00 

•f 

406  1  7 

5.34 

1.19 

70 

0.20 

43  ±2 

130120 

1.51 

1.67 

70 

0.40 

*t 

192120 

2.23 

1.45 

70 

0.70 

tf 

360  1  40 

4.14 

1.24 

70 

1.00 

** 

390 1 70 

4.51 

1.22 

90 

0.20 

60  ±5 

120 1 50 

1.02 

1.98 

90 

0.40 

91 

410 1 30 

3.38 

1.29 

90 

0.70 

•f 

6401  100 

5.28 

1.18 

90 

1.00 

II 

800  1  200 

6.31 

1.16 

rfV 

a.  Number  of  coulombs  of  Au(I)  reduced  over  a  3.3-3.  electrode  area. 


b.  Determined  from  transmission  electron  microscopy. 

c.  Defined  as  the  particle  length,  A.  divided  by  the  particle  diameter,  B.  In  the  text.  A  = 
2a  and  B  =  2b. 

d.  M-G  screening  parameter  (see  text). 
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Table  n. 


Comparison  of  Xmax  Values  from  Experiment,  me  Theory  and  the 
Dynamical  M-G  Model  (Xoas  In  nm)  for  8pherical  Particles 


Particle 

Radius 

(nm)* 

Experimental 

Values1* 

me 

Theory 

Calculation 

Dynamical 

M-G  Model 

Calculation 

30 

556 

556 

568 

38 

559 

576 

608 

43 

604 

592 

636 

60 

656 

* 

658 

756 

a.  Determined  from  transmission  election  microscopy. 

b.  Extrapolated  values  from  experimental  data  shown  In  Figure  6  (see  text). 


Figures 
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Figure  1.  Transmission  electron  micrographs  of  facial  microtome 
sections  of  gold  plated  anodic  alumina  membranes  grown  under  the 
conditions:  (A)  40  V  in  4%  oxalic  acid;  (B)  50  V  in  4%  oxalic  add;  (C) 
70  V  in  2%  oxalic  add;  (D)  90  V  in  2%  oxalic  add.  Scale  bar  in  each 
case  is  20  nm.  Note:  gold  nanopartldes  are  not  visible  in  these  sections. 


Figure  2.  Schematic  of  optical  transmission  experiment  for 
gold/alumina  composites.  The  cylindrical  partides  are  aligned  with  their 
prindpal  axes  parallel  to  the  direction  of  light  incidence. 


Figure  3.  Transmission  electron  micrographs  of  gold  partides  prepared 
by  electrodeposition  of  gold  into  oxalic  add  membranes  prepared  under 
40  V  anodizing  conditions.  The  amount  of  gold  in  terms  of  the  number  of 
coulombs  of  Au(I)  reduced  per  unit  area  of  pores  in  the  membrane  is 
given  with  each  image.  Scale  bar  in  each  case  is  20  nm. 

Figure  4.  UV/Vis/NIR  spectra  for  blank  anodic  alumina  membranes 
and  gold  partide/ alumina  composites  in  which  the  partides  are  nearly 
spherical.  Curve  1:  partide  radius  60  nm;  Curve  2:  partide  radius  43 
nm;  Curve  3:  partide  radius  30  nm;  Curve  4:  blank  30  nm  pore  radius 
alumina;  Curve  B:  blank  43  nm  pore  radius  alumina;  Curve  6:  blank  60 
nm  pore  radius  alumina. 
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Figure  5.  UV/Vis/NIR  spectra  for  gold  particle / alumina  composites 
containing  particles  of  various  aspect  ratios.  Particle  length  in  all  cases 
equals  2a.  (A)  30  nm  pore  radius  systems:  Curve  1.  Blank  membrane 
(ccl  30  pm  in  thickness)  prepared  at  40  V  in  4%  oxalic  add;  Curve  2. 

2a  =  92  nm;  Curve  3.  2a  =  132  nm;  Curve  4.  2a  =  268  nm. 

(B)  60  nm  pore  radius  membranes  prepared  at  90  V  in  2%  oxalic  add: 
Curve  1.  Blank  alumina  membrane;  Curve  2.  2a  =  120  nm;  Curve  3. 
2a  =  410  nm;  Curve  4.  2a  =  800  nm. 


Figure  6.  Summary  of  Xmar  dependence  on  particle  shape.  The 
cylindrical  radii  are  given  next  to  the  curves.  The  values  for  the  aspect 
ratio  (A/B)  are  taken  from  Table  1.  The  M-G  calculation  yields  the  Xmax 
values  shown  in  the  curve  labeled  M-G. 

Figure  7.  Spectra  calculated  from  Mie  Theory  (Equations  1-3).  Input 
parameters  as  follows:  Curve  1:  radius  a  ■  30  nm.  partide  density 
N  =  57.5  pm*1.  Curve  2:  43  nm,  N  =  19.5  pm*1.  Curve  3:  a  =  60  nm, 

N  =  9.95  pm'1.  The  optical  path  length  is  1  pm  in  all  cases.  Host 
medium  no  s  1.6.  . 

Figure  8.  Absorption  spectra  calculated  from  simple  M-G  Theory  using 
Equations  6-8  (Inset:  The  three  layer  model  assumed  for  calculations). 
The  input  parameters  are  as  follows:  Curve  1.  K  =  2,  fm  =  0.0065,  d2  -  1 
mm,  di  ■  28  pm;  Curve  2.  K  =  1.4,  fm  =  0.025,  d2  =  1.07  pm. 
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di  *  27.93;  Curve  3.  K  *  1.22,  fm  =  0.034,  d2  *  1.21  4m.  di  =  27.79  4m. 
The  da  layer  equals  1  tun  in  all  cases. 


Figure  9.  Absorption  spectra  calculated  using  M-G  Theory  with  size- 
adjusted  depolarization  factor  (A)  30  nm  radius  system  input 
parameters:  Curve  1.  a  ■  46  nm.  fm  *  0.0065.  da  *  1  pm,  di  =  28  nm; 
Curve  2.  a  ■  66  nm.  fm  *  0.0133.  da  *  1.070  pm.  di  »  27.93  pm; 

Carve  3.  a  *  0.130,  fm  -  0.217,  da  ■  1.206  pm.  di  *  27.79  pm.  (B)  60 
nm  radius  system  input  parameters:  Carve  1.  a  ■  62  nm.  fm  =  0.0093. 
d2  *  1.0 pm.  di  *  28  pm;  Carve  2.  «> 206 nm.  fm  «  0.024.  d2  =  1.288 
pm,  di  «  26.71  pm;  Carve  3.  a  *  385  nm,  fm  *  0.0314,  d2  *  1.520  pm, 
di  *  26.48  pm.  The  da  layer  equals  1  pm  in  all  cases. 


Figure  10.  Extinction  maxima  calculated  using  M-G  Theoiy  with  size- 
adjusted  depolarization  factor  plotted  against  the  natural  logarithm  of 
the  particle  aspect  ratio  (A/B). 
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